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A simple "trick" is proposed, which allows to perform exactly the site-averaging procedure required 
when developing integral equation theories of interaction site models of macromolecular fluids. It 
shows that no approximation is involved when the number of Ornstein-Zernike equations coupling the 
site-site correlation functions is reduced to one. Its potential practical interest for future theoretical 
developments is illustrated with a rederivation of the so-called molecular closures. 



In the recent years, the application of liquid state in- 
tegral equation methodsi to soft matter systems allowed 
significant advances in our understanding of their prop- 
erties. One of the most prominent developments in this 
field is the formulation of the Polymer Reference Interac- 
tion Site Model (PRISM), extending to macromolecular 

\ systems the Reference Interaction Site Model (RISM) of 
Chandler and Andersen?? Mostly applied to polymeric 
fiuidsr^ it has recently been considered to study fluids of 
platelet-like colloids 

It is usually derived from the RISM theory as follows. 
Consider a one-component fluid with molecules made of 
N identical sites, i.e. the interaction potentials between 
any pair of sites on different molecules are the same. 
These sites are in general non-equivalent (except those 

\ related by obvious symmetries) and, at the pair level, 
the structure of the fluid has thus to be characterized by 

\ a matrix oi N x N site-site total correlation functions 
hij{r), where indices i and j refer to types of sites on 
different molecules. The RISM approach relies on the 
generalized Ornstein-Zernike matrix equation 

! ^ 

I ^iji"^) = ^ Wife * Cki * [uJij + phij]{r), (1) 

I fc,;=i 

where p is the molecular number density and * denotes 
convolution integrals. Cy (r) is the direct correlation func- 
tion between sites of types i and j on different molecules, 
and u!ij{r) the intramolecular distribution function of 
' sites of types i and j. To solve Eq. it has to be 
supplemented by some appropriate closure. In the case 
of hard-sphere sites of diameter d, it can be the Percus- 
Yevick approximation hij{r) = — 1 if r < d, Cij{r) = if 
r > 

PRISM is obtained by making the simplifying assump- 
tion that all N sites are equivalent. This is true for ring 
molecules for instance, but it amounts to a neglect of end 
effects for linear molecules or of edge effects for platelets. 
Then, using hijijr) ~ hir) and Cij{r) ~ c(r) for all j,j's, 
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one finds that Eq. |^ reduces to the single scalar equa- 
tion 

h{r) =uj'i'C^[uj + pNh]{r), (2) 

with 

1 ^ 

= ]^ 51 ujtjir). (3) 

tj(r) is the total intramolecular distribution function. 
Using a perturbation theory approach, Curro and 
Schweizer^ have shown that the best choices for h{r) and 
c(r) are the site averaged values, 

1 ^ 

and similarly for c(r), thus demonstrating that PRISM 
involves a pre-averaging of chain end (or platelet edge) ef- 
fects and evidencing the practical relevance of the theory, 
since the doubly summed quantity of Eq. Q , to which it 
gives access, appears naturally in the calculation of the 
scattering properties or of the thermodynamics of the 
system. 

In the present Note, I want to put forward a "trick" for 
the derivation of PRISM from RISM equations, showing 
that the reduction of the N x N matrix equations of the 
latter to the scalar equation of the former can in fact be 
done without invoking any approximation. 

The basic idea of this derivation is that, if all sites on 
a molecule are identical in term of their interaction po- 
tential, then there is a permutation symmetry between 
them. To take this symmetry into account, one can ar- 
tificially decompose each molecule into two objects: a 
backbone made of N sites, labelled 1 to N, in which are 
embedded the geometry and the intramolecular interac- 
tions of the original molecule, and a set of N identical 
"atoms" , labelled 1 to N, carrying the intermolecular in- 
teractions and occupying at random the backbone sites 
such that all permutations of N atoms have equal proba- 
bility. Clearly, this delocalisation of the interaction sites 
does not correspond to any real physical process and is 
only introduced for formal consideration. 

One can write RISM equations for the mobile atoms, 
which have the usual form now with hatted indices. 
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These simply are formal definitions of the direct correla- 
tion functions (r) , applicable to any molecule within a 
site-site description. 

In this system, all interaction sites are by construction 
equivalent. Thus, for all i and j, one has h^ir) ~ h'[r), 

cij{r) = c'(r), and = ^'('')- addition, 

it is obvious that the total distribution functions of the 
present model coincide with the site-averaged quantities 
obtained within the conventional RISM approach, i.e. 
h'{r) = h{r) and uj'{r) = uj{r). It immediately follows 
from Eq. |^ with hatted indices that Eq. (0, with c re- 
placed by c', is exactly obeyed. Note that nothing can 
be said a priori on the relation between c and c' and 
in particular there is no reason to think they should be 
equal. 

So, what did we learn ? From a practical point of view, 
not much: within the present approach, the reduction of 
RISM to PRISM has been made a tautology and, owing 
to the odd nature of the backbone plus delocalized atoms 
system, the final situation is somewhat equivalent to sim- 
ply taking Eq. JSJ as a formal definition of the total direct 
correlation function. The problem can be considered as 
only being displaced to the question of the non-trivial 
relation between the c^ 's of the usual RISM approach 
and the cj-'s of the present fictitious system. However, 
the fundamental point remains: the reduction from N'^ 
to 1 of the number of Ornstein-Zernike equations does 
not involve any approximation and, at variance with the 
usual claim, nothing like end or edge effects is lost when 
going from the RISM level of description to the simpler 
PRISM one. As usual with integral equation theories, if 
PRISM, compared to a full RISM treatment, can appear 
inadequate in some case, it is only because of our inabil- 
ity to guess a correct expression for the direct correlation 
function. 

More interesting is probably the fact that the present 
picture could provide new insights on the PRISM formal- 
ism and thus lead to improvements of the theory. To il- 
lustrate this point, I propose a derivation of the so-called 
molecular closures^ developed for models where the in- 
termolecular site-site interaction consists of a hard core 
part of diameter d plus a soft tail potential. 



Starting with the classical RISM approach, Cy (r) for 
r > d is often given the interpretation of a renormal- 
ized intermolecular site-site interaction potential. In the 
backbone plus delocalized atoms system, it can thus be 
considered as an effective interaction between backbone 
sites. It is effective because the backbone sites do not 
directly interact and the interaction is actually mediated 
by the delocalized particles. Consider a configuration in 
which particle k on one molecule and particle I on an- 
other molecule occupy respectively the backbone sites i 
and j. Using the simplest superposition approximation, 
the contribution to Cij{r) of two particles, rh on the same 
molecule as k and n on the same as I, can be approx- 
imated as Wj,^ * Cmfi * ijj^ilr). Summing over all mo- 
bile particles and averaging over all their permutations, 
one eventually finds Cij{r) ~ tu * c' * w(r) for r > d& 
Then, using a so-called atomic closure for c'(r), Eq. Q 
with c replaced by c', and the impenetrability constraint 
h{r) = —1 ii r < d, we are left with the same equations as 
those proposed in Ref. 0, but with a somewhat different 
physical content. This points to possible improvements 
of the molecular closures: indeed, it has for instance been 
shown that in some cases, simply convoluting by uj is not 
the best way to give an account of the macromolecular 
connectivity.^ 

Finally, like in PRISM, equivalence of all sites on a 
molecule is often assumed in the framework of the density 
functional theory, in order to derive equations for macro- 
molecular fluids from theories of heteronuclear site-site 
fluids?- - One could thus think of using the present "trick" 
in this context. However, since it is already known 
that the excess free-energy functional can exactly be ex- 
pressed in terms of the total monomer density''^^ (which 
can be thought of as a site-averaged quantity for inho- 
mogeneous systems), here again the proposed approach 
would probably not lead to significant practical progress, 
but only help clarifying the formal relation between the 
site-specific and site-averaged levels of description. 
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